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Abstract of DE4431 349 

A method for adaptive digital image signal 
filtering for fast magnetic resonance (MR) image 
processing is provided. A minimum modification 
to the original image structures is tolerated, while 
noise smoothing is carried out at a prescribed 
filtering level. The filtered image is a weighted 
combination of four sub-images which have been 
obtained from the low-pass filter of the original 
image along four main directions. The weighting 
of each sub-image is controlled by the 
differences owing to these sub-images and the 
original image. The resulting image is then 
focused asymetrically in order to magnify the 
boundaries of the image structure. The overall 
effect of this filter structure is an effective 
adaptive noise reduction and edge magnification 
with efficient execution accompanied by the use 
of matrix computers. The high regulatity and the 
high parallelism of the process render the latter 
suitable for efficient execution in conjunction with 
the use of high grade integrated (VLSI, very large 
scale integrated) circuits or multiprocessor 
systems. The efficiency of the filter for effective 
reduction of image noise and the protection of 
important image structures are discussed and 
demonstrated with the use of a plurality of MR 
images from an MR imaging system having a low 
field strength. 
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@ Adaptives Filtern von digitalen Bildsignalen 

(57) Es wird ein Verfahren zum adaptive n digitalen Bildsignal- 
filtern fur sen n ell e Magnetresonanz (MR)-Bildverarbeitung 
vorgesehen. Eine minimale Modification an den ursprungii- 
chen Bildstrukturen wlrd geduldet, wahrend die Reuschglat- 
tung auf einem vorgegebenen Rfterniveau durchgefuhrt 
wird. Das gefilterte Biid 1st eine gewlchtete Komblnation aus 
vier Unterblldern, die aus dem Tiefpafif ilter des ursprOngli- 
chen Bildes entlang vier Hauptrichtungen erhalten wurde. 
Die Gewichtung an jedern Unterbild wird gesteuert durch die 
Unterschiede durch diese Unterbilder und das ursprungliche 
Bild. Das resultierende Bild wird dann nicht-symmetrisch 
gescharft, urn die Bildstrukturgrenzlinien zu verstarken. Der 
Gesamteffekt dfeser Filterstruktur ist eine effektive adaptive 
Rauschreduzierung und Kantenverstarkung mlt einer effi- 
zienten Durchfuhrung unter Verwendung von Matrizenrech- 
nem. Die hone Regelmafclgkeit und der hohe Prallellsmus 
des Vorgangs machen diesen geeignet fur eine effiziente 
Durchfuhrung unter Verwendung hochgradig integrierter 
{VLSI, very large scale integrated) SchaJtungen oder Multi- 
Prozessorsysteme. Die Leistungsfahlgkeit des Filters fur die 
effektive Reduzierung von Bildrauschen und die Bewahrung 
wichtiger Bildstrukturen wird diskutJert und demonstriert 
unter Verwendung mehrerer MR-Bilder von einem MR-Ab- 
bildungssystem mit geringer Feldstarke. 
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Beschreibung 

Die vorliegende Erfindung betrifft im allgemeinen das adaptive Fiitern von digitalen Bildsignalen, die die 
angeordneten Bildpunkte einer digitalen Abbildung darstellen. Die Erfindung findet insbesondere, aber nicht 
5 ausschlieBlich, Anwendung auf die Rauschglattung und die Strukturscharfung in medizinischen Bildsignalanord- 
nungen von der Art, die von Magnetresonanzabbildungs (MRI)-Systemen erzeugt und dargestellt werden. 

Digitale Bildsignale (welche die angeordneten Bildpunkte einer auf einer Kathodenstrahlrdhre oder einer 
anderen Form von Bildsignal-Transducer darstellbaren digitalisierten Abbildung darstellen) sirid derzeit aus 
vielen Quellen einfach erhaltlich. Auf dem Gebiet der medizinischen Abbildung werden solche Abbildungen 
io typischerweise von CT-Scannern und MRI-Systemen erzeugt Seibst traditionelle medizinische Rontgenabbil- 
dungen, die anfanglich auf photographischem Film gemacht wurden, konnen abgetastet und in ein digitalisiertes 
Bildsignalformat umgewandelt werden. 

Digitale Bildsignale weisen typischerweise eine ordentliche Anordnung von digitalen mehrfach-Bit-Daten- 
werten auf, von denen jeder Wert (z. B. Worter oder "Bytes") die Intensitat (und/oder Farbe) eines bestimmten 
15 vorgegebenen Bildpunktes in einem (z.B. durch Rasterabtastung und Umwandlung der Signale in visuelle 
Abbildungen) anzuzeigenden visuellen Abbildung darstellt 

Leider weisen solche digitale Bildsignale typischerweise Rauschkomponenten auf, welche der Bildqualitat 
sowohl in analytischer als auch asthetischer Hinsicht abtragltch sind. Insbesondere im Fall medizinischer Abbil- 
dungen sind solche Rauschprodukte storend, da sie zu ungenauen medizinischen Diagnosen fQhren oder eine auf 
20 der betrachtbaren Abbildung basierende Diagnose unmoglich machen. Demzufolge wurden betrachtliche An- 
strengungen aufgeboten, urn Rauschkomponenten aus digitalen Bildsignalen herauszufiltern, so daB die resultie- 
rende gefilterte Abbildung bei der Darstellung eine analytisch genauere und asthetisch befriedigeridere Darstel- 
lung des in der Anzeige dargestellten Gewebes eines Patienten darstellen. 

Obwohl viele ziemlich komplizierte Bildfilter vorgeschlagen und/oder angewendet wurde, haben viele von 
25 ihnen eine begrenzte Anwendbarkeit und/oder sind andererseits oft ziemlich unpraktisch anzuwenden in einer 
typischen klinischen Umgebung, in der Zeit und Geratschaften zwangsweise eingeschrankt sind (z. B. im Ver- 
gleich mit der Umgebung der Laborforschung, in der gewissermaBen unbegrenzte Zeit und Geratschaften fur 
wiederholte Versuche der Bildverstarkungsfilterung bereitstehen konnen). Demzufolge besteht immer noch ein 
betrachtliches BedQrfnis nach adaptiven Filterungstechniken fur digitale Bildsignale, die sowohl hochef fektiv als 
30 auch wirtschaf tlich (im Hinbiick auf die erforderliche Filterverarbeitungszeit und die erforderlichen Filtergerat- 
schaften) sind. 

Auf dem Gebiet der Magnetresonanzabbildung (MRI) erhalt man typischerweise eine zweidimensionale 
digitalisierte Abbildung von einem Magnetresonanzabbildungssystem, welches ffir die medizinische Diagnose 
ein Schnittbild von NMR-Nuklei darstellt, welche Gewebestrukturen innerhalb eines dreidimensionalen Objek- 

35 tes definieren. Leider sind die meisten MR- Abbildungen durch verschiedene Verzerrungen verschlechtert 
Manchmal resultieren diese Verzerrungen in ungeeigneten Darstellungen eines Objektgewebes und machen 
eine genaue Bildanalyse schwierig. Unter diesen Verzerrungen ist die Rauschverschmutzung oft die gravierend- 
ste. Effektiye Korrekturen weiterer Verzerrungen stehen oft in engem Zusammenhang mit erfolgreicher Rau- 
schreduzierung und basieren darauf. Daher ist die Rauschreduzierung ein zentrales Thema bei der MR-Bildver- 

40 arbeitung fur bessere Bildqualitat 

Einige Probleme und Erfordernisse bei der MR-Bildverarbeitung, die sich im allgemeinen auf deren Hinter- 
grund beziehen, finden sich in den foigenden Schriften (auf die hiermit vollinhaltlich Bezug genommen wird): 

F.W. Wehrli, Biomedical Magnetic Resonance Imaging- Principles, Methodology, and Applications, Hrsg. 

45 F.W. Wehrli, D. Shaw and J.B. Kneeland, VCH Publishers, Ino, 1988; 

L. Kaufmann, L£. Crooks and J. Carlson, Technology Kequirements for Magnetic resonance Imaging 
System^, in: Proceedings of Technology Requirements for Biomedical Imaging, IEEE Computer Society 
Press, May 1991; V 
AJC Jain, Fundamentals of Digital Image Processing, Prentice Hall, Englewood Cliffs, N J 07632, 1989; 

so WJCPratt, Digital Image Processing, John Wiley & Sons, Ino, 1987; R.C Gonzalez and P. Wintz, Digital 

Image Processing, Addison-Wesley Publishing Company, 1988; 

JS. Lee, "Digital Image Enhancement and Noise Filtering by Use of Local Statistics", IEEE Trans. Patt AnaL 
Machine IntelL, Vol PAMI-2, No. 2, S. 165— 68, Marz 1980; 

J.W. Woods, S. Dravida and R. Mediavilla, "Image Estimation Using Doubly Stochastic Gaussian Random 
55 Field Models", IEEE Trans. Pattern. AnaL Machine Intel!. Vol PAMI-9, S. 245 -53, MSrz 1987; 

L.L. Scharf, Statistical Signal Processing: Detection, Estimation, and Time Series Analysis, Addison-Wesley 
Publishing Company, 1990; 

L.B. Jackson, Digital Filters and Signal Processing, Kluwer Academic Publishing, 1 989. 

60 Das Rauschfiltern von Abbildungen ist im wesentlichen ein GlattungsprozeB, aber einfaches TiefpaBfiltern 
verzerrt typischerweise Kanten bzw. Umrisse und weitere Gewebestrukturen und schadigt damit die Treue der 
Strukturgewebeabbildung. Diese Gewebestrukturen sind fur den menschlichen Betrachter signifikant Sie sind 
auch sehr wichtig fQr die automatische Bildanalyse und -synthese, wie etwa Computer-Sicht- und Bildregistrie- 
rung. Einige Filter wurden vorgeschlagen, urn dieses Problem erfolgreich anzugehen, wie adaptives Fiitern auf 

65 der Grundlage von lokalen Bildstatistiken (z. B. siehe J.S. Lee, "Digital Image Enhancement and Noise Filtering 
by Use of Local Statistics", IEEE Trans. Patt Anal. Machine Intell„ Vol. PAM1-2, No. 2, S. 165-68, Marz 1980), 
und Kalman-Filtertechniken unter Verwendung von Mehrfachbildmodellen (z. B. siehe J.W. Woods, S. Dravida 
and R. Mediavilla, "Image Estimation Using Doubly Stochastic Gaussian Random Field Models", IEEE Trans. 
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Pattern. AnaL Machine Intell. VoL PAM1-9, S. 245—53, Marz 1987). Da fur medizinische Diagnoseabbildung eine 
Studie der MR oft eine Sequenz einer Mehrzahl von Bildern erfordert, und da jedes digitale Bild eine groBe 
Anzahl von angeordneten Bildpunkten enthait, ist die Verarbeitungseffizienz des Filterns eine grundlegende 
ErwSgung bei der Wahl von Filtern fur eine schnelle medizinische Bildverarbeitung. Somit, obwohl letztendiich 
viele Filter eine akzeptable Filterleistung zur VerfQgung stellen mogen, schrankt eine niedrige Verarbeitungsge- 
schwindigkeit oft ihre breite Anwendung bei der medizinischen Abbildung ein. 

Die vorliegende Erfindung stelit ein insbesondere effektives und effizientes Filter fQr die MR- Bildverarbei- 
tung zur VerfQgung. Eine Mehrzahl von Satzen digitalisierter Unterbildsignale wird durch das TiefpaBfiltern 
eines eingegebenen Satzes von digitalen Bildsignalen entlang einer jeweils entsprechenden Mehrzahl von 
Richtungen innerhalb des Bildes erzeugt Danach wird ein Satz von rauschgeglatteten Bildsignalen als eine 
gewichtete Kombination dieser ersten Mehrzahl von Satzen von Unterbildsignalen erzeugt Vorzugsweise ist 
die dem Beitrag jedes Unterbildes zugeordnete Gewichtung eine Funktion des Unterschiedes zwischen dem 
Unterbild und dem ursprQngiichen Bild, wodurch das resultierende Filter adaptiv fQr die Gewebestruktur 
gemachtwird 

Bei dem bevorzugten beispielhaften Ausfuhrungsbeispiel werden die Kanten von Gewebestrukturen inner- 
halb des rauschgeglatteten digitalen Bildes gescharft, indem anschlieBend aus dem Satz von rauschgeglatteten 
digitalen Bildsignalen eine Mehrzahl von weiteren digitalisierten Satzen von Unterbildsignalen durch HochpaB- 
filtern einer jeweils entsprechenden Mehrzahl von Richtungen (die die gleichen wie die fur das vorhergehende 
TiefpaBfiltern sein konnen oder auch nicht) erzeugt wird, und anschlieBend ein ausgegebener Satz von digitalen 
Bildsignalen als eine gewichtete Kombination einer solchen Mehrzahl von weiteren Satzen von Unterbildsigna- 
len erzeugt wird. Dieses Kantenscharfungsfilter kann auch verwirklicht werden durch Verwendung der oben 
erwahnten Unterschiede zwischen tiefpaBgefilterten Unterbildern und dem ursprunglichen Bild — womit erfor- 
derliche Filterressourcen weiter erhalten bieiben. 

Bei einem beispielhaften AusfQhrungsbeispiel wird sowohl das TiefpaB- als auch das HochpaBfiltern unter 
Verwendung von LSI (linear space invariant)-Filtern verwirklicht Dennoch werden strukturadaptive Filteref- 
fekte erzeugt durch darauffolgendes Erzeugen einer adaptiv gewichteten Kombination solcher LSI-Filteraus- 
gange. Eine solche Verwendung vori LSI-Filtergruppen zur Verwirklichimg von strukturadaptiver Bildverarbei- 
tung stelit nicht nur eine effektive Rauschunterdrttckung, sondern auch eine Strukturverbesserung von MR- Bil- 
dern zur VerfQgung, wahrend sie gleichzeitig effektive Anwendungen unter Verwendung eines Matrixprozes- 
sors, paralleler Signalverarbeitungskanale und dergleichen erlaubt 

Obwohl die LSI-Filter nahezu in jeder beliebigen Anzahl von unterschiedlichen Richtungen betrieben werden 
konnen, die die GroBe der digitalen Bildanordnung erlaubt, wendet eine gegenwartig bevorzugte effektive und 
effiziente Durchf Qhrung das LSI-Filtern nur in vier Richtungen an, die in Winkelintervallen von im wesentlichen 
45° innerhalb des Bildraumes ausgerichtet sind Ein Beispiel eines geeigneten LSI-TiefpaBfilters erzeugt fur 
jeden Bildpunkt eine gewichtete Kombination einer vorbestimmten Anzahl seiner nachsten Nachbarn entlang 
einer Linie, die parallel zu der Richtung dieses bestimmten Filters ist Eine beispielhafte LSI-HochpaBfilterung 
ist vorzugsweise nicht-symmetrisch im Hinblick auf sowohl die raumliche Richtung als auch die Polarit&t — und 
kann auch bewirkt werden als eine gewichtete Kombination einer vorbestimmten Anzahl von nachsten Nach- 
barn entlang einer Linie parallel zu der Filterrichtung. 

Das beispielhafte Filter stelit eine auBergewohnlich schnelle Magnetresonanz (MR)-Bildverarbeitung zur 
VerfQgung, indem es bewirkt, was wahrscheinlich so gut wie eine minimale Modifikation an den ursprQngiichen 
Bildstrukturen ist, wahrend es jedoch eine Rauschglattung auf einem vorgegebenen Filterpegel durchfQhrt Das 
gefilterte Bild ist eine gewichtete Kombination der Bildpunkte in vier Unterbildern, die durch das TiefpaBfiltern 
des ursprlinglichen Bildes entlang vier Hauptrichtungen erhalten sind Die Gewichtung an jedem Unter bild-Bild- 
punkt ist gesteuert durch die Unterschiede zwischen bildpunkten in diesen Unterbildern und einem entsprechen- 
den Bildpunkt im ursprOngltchen Bild Das resultierende Bild ist dann nicht-symmetrisch gesch&rft, urn die 
strukturellen Grenzlinien von abgebildetem Gewebe zu verstSrken. Der Gesamteffekt stelit eine adaptive 
RauschunterdrOckung und eine Kantenverstarkung mit einer effizienten Durchf uhrung unter Verwendung von 
Matrixprozessoren zur VerfQgung. Die hohe Regularitat und der hohe Parallelismus des Filters macht ihn auch 
geeignet fQr eine DurchfQhrung unter Verwendung von VLSI (very large scale integrated)-Schaltungen oder 
Multi-Prozessoren-SystemeiL Die Leistung des Filters bei der UnterdrQckung von Bildrauschen und der Bewah- 
rung von wichtigen abgebildeten Gewebestrukturen ist weiter unten unter Verwendung mehrerer MR-Bilder 
von einem MR- A bbildungssys tern mit geringer Felds tarke demonstriert 

Wie viele weitere adaptive Filter versucht das beispielhafte Filter dieser Erfindung, Bildrauschen zu unter- 
drucken und doch wichtige Bildstrukturen zu erhalten. Hier wird dieses Ziel der Minimierung einer Modifikation 
am ursprUnglichen echten Bild (d h. dem theoretisch rauschlosen, perf ekten Bild) angegangen durch die Auswahl 
von Fllterparametern, die ein vorbestimmtes eingeschranktes Glattungsniveau in homogenen Bereichen von 
Bildern erzeugen. Die bevorzugten AusfUhrungsbeispiele vermeiden eine herkommiiche Bildanalyse, wie sie 
typischerweise bei adaptivem Filtern verwendet wird (z. B. wie etwa Kantendetektion unter Verwendung eines 
Gradient en oder Laplace-Operators, und Bildmodellidentifikation) und welche resourceverzehrende Berech- 
nungen umfassea Das anfangliche rauschgeglattete, gefdterte Bild wird hergestellt als eine gewichtete Kombi- 
nation von vier Unterbildern, die durch TiefpaBfiltern des ursprunglichen Bildes entlang vier Hauptrichtungen 
erhalten wurden. Diese TiefpaBfiltervorgange in jeder Richtung sind nicht-adaptiv und konnen daher sehr 
effektiv unter Verwendung von Matrixprozessoren durchgefQhrt werden. Das resultierende Bild stelit jedoch 
eine adaptive Glattung mit der erforderlichen RQcksicht auf die Bewahrung echter Gewebebildstrukturen zur 
VerfQgung, da Gewichtungsparameter fQr die Synthetisierung des endgQltigen TiefpaBbildes durch die Unter- 
schiede zwischen diesen Unterbildern und dem ursprunglichen Bild gesteuert sind Diese Unterschiede stellen 
Informationen Qber Bildstrukturen zur VerfQgung und stellen damit eine strukturadaptive Filterung zur VerfQ- 
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gung. Daher konnen Bildstrukturen erhalten werden, wahrend das Bildrauschen effektiv verririgert wird Des ' 
weiteren, da der Filtervorgang einen stark parailelen Aufbau aufweist, kann er unter Verwendung mehrerer 
synchronisierter Mikroprozessoren oder VLSI-Schaltungen durchgefQhrt werden, urn die Geschwindigkeit der 
Verarbeitung weiter zu erhohen. Als Ergebnis kann das Filter f ilr eine schnelle Rauschunterdruckung bei MR 
s und weiteren medizinischen Abbildungsanwendungen niitzlich sein. Vorzugsweise wird auch ein nicht-symme- 
trisches Kantenscharfungsfilter vor dem Vorsehen eines endgQltigen zusammengesetzten Filterausgangs fur die 
Bildstruktur-Grenziinienverstarkung verwendet 

Diese wie auch weitere Aufgaben und Vorteiie der vorliegenden Erfindung sind vollstandiger zu verstehen 
uhd zu wurdigen anhand der folgenden, detaillierteren Beschreibung gegenwartig bevorzugter AusfQhrungsbei- 
10 spiele der vorliegenden Erf indung in Verbindung mit der beigefugten Zeichnung. 

Eszeigt 

Fig. 1 einen schematischen Oberblick Qber ein MRI-System, das derart modifiziert ist, daB es Filtern gemaB 
einem Beispiel der vorliegenden Erfindung zur VerfQgung stellt; 

t 2 veraaischaulicht eine beispielhafte Ausrichtung von vier Hauptfilterrichtungen zur Verwendung mit 
is LSI^HochpaB. und TiefpaBfiltern gemaB dem Beispiel von Fig. 1 der vorliegenden Erfindung; 

Fig. 3 eine schematische Abbildung eines beispielhaf ten adaptiven TiefpaB-Bildfilters zur Verwendung mit der 
vorliegenden Erfindung; 

Fig. 4 eine schematische Abbildung eines nicht-symmetrischen KantenschSrfungsfilters zur Verwendunc in 
der vorliegenden Erf indung; 

20 ^ Fig. 5 ein vereinfachtes Ablaufdiagramm einer geeigneten Computerprogramm-Subroutine, die in dem MRI- 
System von Fig. verwendet werden kann, um ein. AusfQhrungsbeispiel des Filterns gemaB der vorliegenden 
Erfindung praktisch durchzufQhren; 
Fig. 6 ein Diagramm eines durchgef Qhrten AusfQhrungsbeispiels der Erfindung; 

Fig. 7 eine Photographie einer eingegebenen sagittalen Kopfabbildung, die mit einem herkommlichen dreidi- 
25 mensionalen MRI-System hergestellt wurde; 

Fig. 8A, 8B, 8C und 8C Photographien zur Veranschaulichung beispielhafter Unterbilder, die aus dem ausge- 
richteten TiefpaB-LSI-FiItern der Kopfabbildung von Fig. 7 entlang der Richtungen 0°, 90° 45° bzw 135° 
resultieren; ' 

Fig. 9 eine Photographie eines rauschgeglatteten Filterbildes, das als gewichtete ^Combination der Unterbil- 

30 der in Fig. 8A—8D erzeugt wurde; und 

Fig. 10A, 10B, 10C und 10D Photographien zur Veranschaulichung zweier Paare von ursprQnglichen und 
getuterten traiwaxialen Kopfabbildungen von einem zweidimensionalen MIR unter Verwendung von ersten 
10DX * 10A ' 10B ^ ^ UntCr Verwendung von zweiten Spinecho-MR-Bilddaten (Fig. 10C, 

35 Pig. 1 stellt ein typisches herkdmmliches MRI-System dar, das so angepaBt wurde, daB es ein AusfGhrungsbei- 
l P ™loTM rhege ^ ^ J^flpdung praktisch durchfQhit Ein Beispiel fur solch ein System ist das MRI-System 
yon Toshiba. Es kann beispielsweise eine ziemlich groBe polarisierende NMR-Magnetstruktur 10 
aufweisen, die cm im wesentlichen gleichformig homogenes polarisierendes NMR-Magnetfeld B 0 innerhalb 
eines PaUentenabbildungsvoIumens 12 erzeugt Ein geeigneter Schlitten 14 wird verwendet, um den gewunsch- 

40 ten Abschnitt der Anatomie des Patienten 16 in das Biidvolumen 12 einzufQhren. Die magnetischen NMR-Gra- 
dienten in B 0 konnen selektiv erzeugt werden durch elektromagnetische Gradientenspulen, NMR-HF-Muta- 
ttonsimpulse konnen in das Gewebe des Patienten im Abbildungsvolumen Qbertragen werden, und NMR-HF- 
Reaktionen konnen yon dem Gewebe des Patienten Qber geeignete HF-Spulenstrukturen empfangen werden, 
wie der .Fachmann wissen wird. Eine bestimmte MIR-Datenerfassungssequenz von solchen magnetischen Gra- 

45 dientemmpulsen, HF-Mutationsimpulsen und NMR-HF-Reaktionen wird herkommlicherweise erhalten mit 
einem MRI-Sequenzkontrulier 18, der die Qbliche Anordnungen von Gfauieniettir eibern 20, HF-Senderschal- 
uS g o n f 2 " nd HF ' Em P fan S erschaltu ngen 24 steuert, die alle auf geeignete Weise mit elektromagnetischen 
HF-Spulen innerhalb des MRI-System-PortalgerQsts verbunden sind. Die empfangenen NMR-HF-Reaktionen 
werden digitalisiert und an einen MRI-Bildprozessor 26 weitergeleitet, der typischerweise einen Matrizen-Rech- 

50 r j. \ geeignete Computerprogramm-Speichermedien 30 (z. B. RAM in Silizium- oder magnetischen 
Medien) aufweist, in denen Programme gespeichert und selektiv angewendet sind, um die Verarbeitung erfaBter 
MR-Daten zu steuern, um digitalisierte Bildanzeigen auf dem Kathodenstrahlrdhren-Terminal 32 zu erzeugen. 
Das Steuertermmal 32 kann auch geeignete Tastenschalter und dergleichen aufweisen, um eine Bedienersteue- 
rung Qber den MRI-Sequenzkontroller 18 und den dazwischengeschalteten unterstutzenden MR-Bildprozessor 

55 26auszu0ben. K 

In Verbindung mit dem Qblichen MRI-Prozessor 26 wird einem Bediender typischerweise ein AuswahlmenQ 
prasentiert Bei dem AusfQhrungsbeispiel der vorliegenden Erfindung ware eine dem Bediener zur VerfQgung 
stehende Wahl die "Verfeinerung" oder Filterung eines MR-Bildes. Bei dem bevorzugten AusfQhrungsbeispiel ist 
em Standardsatz von Filterparametern fQr den allgemeinen Gebrauch vorgesehen. Falls gewQnscht, kdnnen 

60 jedoch andere Satze von Filterparametern gewahlt werden (moglicherweise aus einer vorhefbestimmten An- 
ordnung von SStzen von Filterparametern), um dem Bediener einige Filteraltemativen zur VerfQgung zu stellen, 
die im Falle einer besonderen Art von medizinischer Abbildung bessere Resultate zur VerfQgung stellen kdnnen. 

Um das Bild durch die Reduzierung von zuftlligen Hochfrequenzrauschkomponenten in dem Bild zu verf ei- 
nem, sind TiefpaBfilter notig. Aber LSI (linear space invariant)-TiefpaBfilter konnen wichtige Bildstrukturen, wie 

65 etwa in Querrichtung ausgerichtete Kanten und Linien, verschwimmen lassen und somit den Bildkontrast 
herabsetzen. Um dieses Problem zu Qberwinden, sind raumlich variante (d h. strukturadaptive) Filter vonndten. 
Im Idealfall sollten solche Filter ein starkes TiefpaBfiltern nur entlang von (d h. parallel zu) Kanten und Linien 
und in homogenen (d. h. nicht-strukturierten) Bereichen von Bildern durchfQhren. Bei dem AusfQhrungsbeispiel 
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wird ein schneller Rauschunterdriickungsvorgang durchgefuhrt durch die Verwendung einer Bank von vier 
ausgerichteten LSI-TiefpaBfiltern. Insgesamt wird dennoch eine strukturadaptive Filterung im endgtiltigen 
Filterausgang verwirklicht durch die Bildung einer gewichteten Kombination der vier Unterbilder von dieser 
Fiherbank. 

Im Fall des ursprUnglichen Bildes xo(ni, n 2 ) von Fig. 2 werden vier TiefpaBfllter verwendet, una vier tiefpafige- 
filterte Unterbilder xi_4 entlang von vier Hauptrichtungen (0°, 45°, 90°, 135°) zu erhalten: 

Gleichungen (1) 

xn(ni,n 2 ) — fii(x<>(ni,n 2 )) 
xi2(n t ,n 2 ) = ffc(xo(nt,n 2 )) 
xi3(ni^i 2 ) = fb(xo(ni f n 2 )) 
xu(ni,n 2 ) «= fi4(xo(ni t n 2 )) 

Die Funktion fu stellt eine TiefpaBfilterf unktion dar, die entlang der i-ten Richtung durchgefuhrt werden soil, 
wenn der Index i «= 1, 2, 3 oder 0, 45, 90 bzw. 135 Grad entspricht Diese Filter konnen als konventionelle 
LSI-TiefpaBfilter verwirklicht werden. Jedes der vier Filter soil die gleiche Filter- w Starke w haben (d. h. wie von 
der Anzahl von Bildpunkten definiert ist, die zur Glattung verwendet wird oder die von der TiefpaBfilter-Band- 
breite spezifiert ist). Wenn beispielsweise vier nachste Nachbarbildpunkte um die Position (ni, n 2 ) zusammen mit 
dem Bildpunkt an dieser Stelle zur TiefpaBfilterung verwendet werden (d. h. insgesamt funf Bildpunkte im 
Filterfenster verwendet werden), kann eine Gruppe von einfachen HefpaBfUtern im Bildpunkt domain durchge- 
fuhrt werden, um vier entsprechende Unterbilder herzustellen: 

Gleichungen (2) 

Die Einstellung der Werte fur die Gewichtungsfaktoren cti andert die Charakteristiken dieser Filter. Fur die 
Linearphase und Normalitat der Filter werden die Gewichtungsfaktoren ai eingeschrankt, so daB sie die 
folgenden Beziehungen aufweisen: 

Gleichungen (3) 



a^j^a., and = 1 

Solche TiefpaBfllter konnen auch im Frequenzbereich unter Verwendung von Matrixprozessoren mit FFT 
(Fast Fourier Transform)-Unterprogrammen durchgefuhrt werden. In dem Fall kdnnen die Frequenzcharakteri- 
stiken durch direkte Spezifizierung von unterschiedlichen Frequenzgangen direkt im Frequenzbereich einge- 
stellt werden. 

Die Unterschiede dieser Unterbilder zum ursprQnglichen Bild werden berechnet als: 
Gleichungen (4) 

di(ni,n2)«=xo(ni,n 2 )— xu(ni,n 2 ) 
d^nt^—xoCnj^)— xi 2 (nun 2 ) 
d 3 (ni^ 2 ) « xo(ni^ 2 )— xi3(ni,n 2 ) 
d4(ni4i2)-=xo(niji 2 )~xi4(niji 2 ) 

Die Absolutwerte dieser Bildunterschiede sind hierbei ausgedruckt als: 

Gleichungen (5) 

33«|d3|,a*4= 3 |d4| 

Das Argument (ni, n 2 ) wurde in diesen Gleichungen weggelassen, um Raum zu sparen, und dies wird aus dem 
gleichen Grund auch in den folgenden AusdrOcken getan. 

Bei Oberlegung ergibt es sich, daB die Werte dieser Bildunterschiede von den Bildgewebestrukturen abhan- 
gen. Wenn TiefpaBfiltern fu(xo(ni, n 2 )) endang von Kanten oder Linien durchgefuhrt wird, werden diese Bild- 
strukturen geglattet und im Ausgang xn(ni, n 2 ) gut erhalten, und daher zeigt der resultierende gefilterte Bild- 
punktwert xu(ni, n 2 ) des Bildes nur einen geringen Unterschied zu dem entsprechenden Bildpunktwert xo(ni, n 2 ) 
des ursprQnglichen Bildes. Ausgange von solchen ausgerichteten TiefpaBfiltern stellen daher eine gut geglattete 
Annaherung an das ursprtagliche Bild in diesen Bildbereichen zur Verfflgung. Wenn jedoch TiefpaBfllter Qber 
Kanten oder Linien hinweg durchgefuhrt wird, dann werden diese Strukturen ernsthaft verzerrt, und daher zeigt 
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der resultierende Bildpunktwert xii(n,, n 2 ) des BUdes einen relativ groBen Unterschied zu dem Bildpunktwert ~ 
Mni, n 2 ) des ursprQnglichen Bildes. Die Werte dieser Unterschiede hangen auch von der "St§rke" (dh der 
Bandbreite) des TiefpaBfilters ab. Starkeres TiefpaBfilter (d. b, engere Bandbreite, groBere Anzahl von Bildpunk- 
ten in emem durchschmttlichen Fenster) resultiert in starker gegtatteten Bildern und damit in grOBeren Unter- 

5 scniedswerten. Da diese Bildunterschiede Informationen aber Strukturen im ursprQnglichen Bild (d. h. Informa- 
tionen uber Annlichkeiten der Unterbilder tnit dem ursprQnglichen Bild) zur VerfQgung stellen, kann man eine 
gegKLttete Annaherung des ursprQnglichen Bildes aus diesen vier Unterbildern synthetisiereri, gesteuert durch 
die Bildunterschiede und eingeschrankt durch einen erwflnschten Glattegrad des gefilterten Ausgangs in homo- 
genen Bildbereichen. Die Glatte eines gefilterten Bildes eines bekannten homogenen Phantoms kann beispiels- 

10 weise bewertet werden durch die Messung von Standardabweichungen (SD, standard deviations) des Bildes in 
diesen bekannten homogenen Bildbereichen. Dieser geforderte Glattegrad bestimmt dann die erforderliche 
mimmale Filterstarke (d. h. die maximale Bandbreite) des TiefpaBfilters. 

Ein strukturadaptiver tiefpaByerarbeiteter Bildausgang xio(n,, n 2 ) kann nun synthetisiert werden als gewichte- 
te Kombination der oben disfcutierten vier Unterbilder: 

15 

Gleichung(6) 

X|0 «= Wll • Xll+Wfe • X12 + W13 • X13 + WI4 • X|4 

20 wobei wii, w a w B . ww (0 <S wr £ 1) Gewichtungsparameter sind, die jeweils den vier tiefpaBgefilterten 
Unterbildern zugeordnet sind, und konnen als Funktionen gii(d,, dz, d 3 , d») der Bildunterschiede bestimmt 
werden Urn die gleichen Durchschwttswerte des tiefpaBgefilterten Bildes xk (n,, n 2 ) und des ursprQnglichen 
Bildes xo(n,, n 2 ) zu behalten, sollten diese Parameter das folgende Verhiltnis aufweisen: 

25 Gleichung(7) 

wntni^+w^ni^+wi^ni^+wi^ni,^) = l,forall(ni,n 2 ) 

1 1 ^ a< ?-, j CT Y oranstehend .en Diskussion uber das Verhaltnis von Bildstrukturen und Bildunterschieden sollte ein 
30 Unterbild nut emem genngeren Unterschied zu dem ursprQnglichen Bild, urn erwunschte Bildstrukturen zu 
erhalten, einen relativ groBeren Beitrag zu dem endgOltigen synthetisierten Bild leisten. Daher sollten die 
Oewicntungsparameter wu(n,, n 2 ) monotone Verringerungsfunktionen des entsprechenden Bildunterschiedes 
di^ni, n 2 ) sein. Erne Form davon ist 
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wobei si ein Skaherfaktor ist, der gewahlt wurde, urn die Gleichung (7) zu erfullen, und wo das hochgesetzte k 
einen exponentionellen Parameter fur das Filter darstellt FQr eine einfache Berechnung von 4* kann der 
Parameter k als positive ganze Zahl gewahlt werden. In Anbetracht der Einschrankung der Gleichung (7) 
konnen die Gewichtungsparameter aus den Bildunterschieden bestimmt werden als 

Gleichung (8) 



= i j-. for i = 1, 2, 3, 4 

IT IT IT IT 

Daher haben diese Gewichtungsparameter die folgenden Beziehungen 
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Gleichungen(9) 
_1_ 

i 

_i_ 

1 

5T 



Durch Einsetzen der Gleichung (8) in die Gleichung (6) kann das synthetisierte strukturadaptive tiefDaBeefil- 
terte Bild eindeutig ausgedruckt werden als 



Gleichung(lO) 



Das Filter bewirkt daher eine reduzierte oder minimale Modifikation an den Daten des ursprfinglichen Bildes 
(und erhalt somit die Bildstrukturen) durch die "AuswahT gewichteter Daten aus den Ausgangen von vier 
ausgenchteten TiefpaBfiltern, wahrend es ein gewisses, vorherbestimmtes Friterungs(Glattungs-)niveau in ho- 
mogenen Bildbereichen beibehalt Wenn beispielsweise der Bildpunkt Xo(n u n 2 ) auf einer Kante oder einer Linie 
entlang der 0° -Richtung ist, dann hat der Ausgang xn(m, n 2 ) aus dem ersten FUter fii(ni, n 2 ) an der Position (n lf n 2 ) 
einen viel kleineren Unterschied zu xtfm, n 2 ) als die Ausgange aiis den drei anderen Filtern. Das heiBt, &i < 3, fur 
\>Z % 3 ; 4 * A / US de f 5 31e . ichun S ( i0 ) er ^t sich an dieser Position x t0 (n,, n2 ) « xu(n,, n 2 ). Gleichzeitig, wenn der 
tfiidpunkt xo(ni, n 2 ) in einem homogenen Bildbereich ist, dann 3j « 3 2 » d 3 « <U, und aus der Gleichung (10) ist 
2? r r en 5f? ,tige Ausgan ^ Xl ° un g efah r %(xii +X12 + X13 + X14) oder der Durchschnitt der Ausgange von alien vier 
TiefpaBfiltern. Somit, obwohl jedes TiefpaBfilter fQr sich nicht adaptiv und jedes ausgegebene Unterbild daher 
entlang seiner jeweihgen Richtungen stark gegiattet ist, ist das Gesamtfilter adaptiv an Bildstrukturen, und das 
synthetisierte Bild ist nur stark gegiattet in homogenen Bereichen und entlang von Kanten und Linien der 
Bildstrukturen. Daher kann das Filter Rauschen effektiv unterdrucken, wahrend es jedoch auch wichtige BUd- 
strukturen erhalt Die Filterleistung kann eingestellt werden durch die Parameter as des TiefpaBfilters ft und k 
des Gewichtungsfaktors wu. Ein groBerer Parameter k macht den Filtervorgang empfindlicher fiir Bildunter- 
schiede und damit empfindlicher fur Bildstrukturen. Starkeres TiefpaBfiltern (d. h. eine solche Wahl von a if daB er 
mehr Bildpunkte zusammen mittelt) in jeder Richtung stellt glattere endgQitige Bilder zur Verfugung. 

Die Form des beispielhaften Filters ist sehr ahnlich zu derjenigen eines Berechners fur die groBte Wahrschein- 
lichkeit ernes Wertes x aus den vier Messungen xii fQr i « 1,23,4, wenn x\ jev/eils Normal -verteilt mit dem 
Mittelwert und den Standardabweichungen 3i sind (z. B. siehe L.L. Scharf, Statistical Signal Processing: Detec- 
t ^ n, . E ?V matlon ' and Time Series Ana,vsis » Addison-Wesley Publishing Co., 1990). Es ist auch interessant, die 
Ahnlichkeit des Bildfiltervorgangs mit Bildinterpolation festzustellen. Man nehme das urspriingliche Bild x(m 
n 2 ) und die vier Unterbilder xu(nj, n 2 ) als fOnf Punkte in einem Bildraum an. Die Distanzen vom Punkt des 
ursprungiichen Bildes zu den vier Unterbildpunkten sind 3i, 3* 3 3 bzw. Aus der Gleichung (10) embt sich, 
wenn k - 1 ist, das endgultige tief paBgefilterte Bild als: 

Gleichung (11) 

111 1 



15 



11 I 1 2 ° 

d% d* d* d* 



m — i.i i — i 



Wie zu sehen ist, ist es in diesem Fall eine lineare Interpolation von vier Unterbildern xii(n u n 2 ), i = 1, 2, 3, 4 
Dies ist die glatteste Annaherung an das urspriingliche Bild aus dem Filtervorgang. Wenn hingegen k — + a: 

Gleichung (12) 

xio(ni,n 2 ) = xu(ni4i 2 X if ai(n lf n 2 )<aj(ni4i2) 
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In diesem Fall ist das tiefpaBgehlterte Bild bestimmt durch einen Interpolationsvorgang des nachsten Nach- 
barn und das Filtern nahert sich einem binaren Selektionsvorgang aus Ausgangen der Filterbank an. Es soli 
angemerkt sein, daB das Filter nun ein direktes Durchschnittsbildungsfilter (wie beispielsweise in A.K. Jain, 
Fundamentals of Digital Image Processing, Prentice Hall, Englewood Cliffs, NJ 07632, 1989 diskutiert ist) fur vier 
Richtungen wird/Obwohl das durch die Gleichung (12) defmierte Filter einen minimalen Bildunterschied ergibt, 
indem es nur ernes der vier Unterbilder wahlt, ist die resultierende Glattung oft nicht adaquat in nicht-texturier- 
ten Bildbereichen. Da nur ein Filterausgang mit dem Minimum 8i fur dieses TiefpaBfiltern verwendet wird, 
werden nur funf (anstatt siebzehn) Bildpunkte verwendet, urn die Rauschglattung durchzufuhren <sogar in 
homogenen Bildbereichen). Eine Erhohung der in der Filtergleichung (12) verwendeten Anzahl von Bildpunkten, 
urn die gleiche Filters tarke wie die von der Gleichung (10) zu erzielen, in nicht-texturierten Bildbereichen, 
beeintrachtigt die Empfindlichkeit des Filters far Bildstrukturen. Somit stellt die Filterstruktur der Gleichung 
(10) erne L6sung des Rauschunterdruckungsproblems mit besserer Erhaltung von Kanten zur Verfugung als das 
durch die Gleichung (12) definierte Filter. * 

Ein Beispiel fur dieses neue FUterschema ist bUdlich in Fig. 2 dargesteilt, wo *u = wr, xu der gewichtete 
is Ausgang von jedem TiefpaBfilterkanal ist Da jeder TiefpaBfilter linear space invariant ist, wird die gleiche 
Operation an vielen Anordnungen von Daten durchgefuhrt Die Operationen zur Bildunterschiedberechnung, 
Gewichtungsparameterbestimmung und endgultige BUdsynthesen sind auch globaler Natur. Daher konnen 
Matnzen-Rechner fur eine effiziente DurchfQhrung dieses Filters verwendet werden. Es zeigt sich, daB dieses 
Filter erne regelmaBige Berechnungsstruktur besitzt und somit fQr VLSI-Hardwaredurchfuhrungen geeignet ist 
20 Das Filter besitzt auch hochgradigen Parallelismus und ist daher in einer Umgebung mit mehrfachen MikrdDro- 
zessorenleicht durchzufuhren. . • 

Wenn Bilder fur den menschlichen Betrachter verarbeitet werden, ist es oft wunschenswert, Bildkanten fur 
asthetisch ansprechendere Ansichten zu scharfen. Der am gewohnlichsten verwendete Ansatz wird als "Un- 
scharfe-"maskierung bezeichnet Es wird durchgefuhrt, indem ein Bild durch einen HochpaBfilter geftthrt wird 
25 und dann ein Bruchteil seines Ausgangs zu dem ursprunglichen BUd wird. Obwohl dieser herkOmmliche Ansatz 
m vielen Fallen ziemhch gute Resultate ergibt, wird auch unerwOnschtes neues Bildrauschen miterzeugt, beson- 
ders wenn das auseegebene Bild einen niedrigen Rauschabstand (SNR) hat Ein weiteres Problem mit diesem 
Ansatz ist es, daB Uberschwingen von hellen Bildkanten hergestellt werden konnen, wenn eine relativ starke 
Scharf ung angewendet wird. 

30 t,3 1t \* nic ^ s ^ met risches Kantenscharfungsverfahren ist daher fur die Verwendung mit der vorliegenden 
Erfmdung bevorzugt, urn diesem Problem entgegenzuwirken. Das Scharfungsfilter ist vorzugsweise nicht-sym- 
metnsch im Hmbhck auf die raumliche Richtung als auch auf das Vorzeichen des Verstarkungssignals. Ein 
Beispiel ra^em solches Kantenscharfungsschema ist in Fig. 4 abgebildet Die Bilddaten xio(ni, ri 2 ) void adaptiven 
TiefpaBfilter von Fig. 3 werden hochpaBgefiltert endang vier Richtungen (die die gleichen seih konnen wie die 
im TlefpaBfiltervorgang des Standes der Technik von Fig. 3, der im wesentlichen unabhangig istl um weitere 
Unterbildausgange xhi(ni, n 2 ), Xh2(n t , n 2 % x h 3(ni, n 2 ) und x h 4(n t , n 2 ) zu erhalten. Das heiBt: 

Gleichungen(13) 

40 Xhi(n lt n 2 ) - fhi(xio(ni4i 2 )) 
Xh2(ni,n 2 ) - fh2(xio(ni ( n 2 )) 
Xh3(ni,n 2 ) — fhs(xio(ni,n 2 )) 
Xh4(ni,n 2 ) » fM(xio(ni^ 2 )) 

45 wo hh fh2, fh3 und fb4 LSI-HochpaBfilter sind, die entlang der Hauptrichtungen von 0, 45, 90 bzw. 135 Grad 
durchgefuhrt werden. Ahnlich zu dem Fall des TiefpaBFilterns k6nnen diese Filter unter Verwendung von 
Matnxprozessoren direkt im Bildbereich oder im Frequenzbereich durchgefuhrt werden. Dann erhalt man die 
gewichtete Summierung der resultierenden vier HochpaB-Unterbilder als 

50 Gleichung (14) 

XhO = Whl • Xhl +Wh2 • Xh2 + Wh3 * Xh3 + Wh4 • Xh4 

Die Gewichtungsparameter ww werden als Funktionen von Ausgangen dieser HochpaBfilter bestimmt und so 
55 emgeschranktdaBihrVerhaltnisist: 

Gleichung (15) 

Whi(ni,n 2 )+wh2(ni,n 2 )-f Wh3(ni^i 2 ) + wh4(ni4i2) = l,for all(m,n 2 ) 

Um Komponenten mit hohen Amplituden aus der kantenQberquerenden Richtung (anders als im Fall des 
TiefpaBfilterns) zu verstarken, sollten diese Gewichtungsparameter einige monoton ansteigende Funktionen der 
Amplituden der entsprechenden UnterbUder sein. Eine einfache Form kann Whi = sh-lxmP sein, wo sn ein 
Skalierfaktor ist, der so gewahlt ist, daB er die Einschrankung von Gleichung (15) erfQilt, und das hochgestellte j 
65 bezeichnet einen exponentiellen Parameter des HochpaBfilters. Der Parameter j kann wiederum fQr einfache 
Berechnungen eine positive ganze Zahl sein. Unter Verwendung dieser Funktion und unter RQcksichtnahme auf 
Gleichung (15) konnen die Gewichtungsparameter bestimmt werden durch 
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Gleichung (16) 

"VT g u {*» *„. x„, x„) = for f . 1b a. 3, 4 

wobei der exponentielle Parameter j verwendet werden kann, urn die Leistung des Filters einzustellea Unter 
Verwendung der Gleichung (16) und der Gleichung (14) lassen sich die kombinierten Daten des hochpaBgefilter- 
ten Bildes eindeutig ausdriicken als: 

Gleichung (17) 



Der Ausgang ww • xhi von jedem Filterkanal ist durch Xhi in Fig. 4 dargestellt 

Das Filter verstarkt selektiv Komponenten mit hoher Amplitude und hoher Frequenz von den HochpaBfiltern 
in Richtungen, die Qber Bildstrukturkanten hinweg (d. h. in Querrichtung dazu) ausgerichtet sind. Filterausgange 
mit hoheren Amplituden stellen groBere Beitrage zum endgultigen Verstarkungssighal zur Verfiigung, und ein 
groBerer Parameter j macht diesen Auswahlvorgang empfindlicher fur Amplituden der HochpaB filterausgange 
in verschiedenen Richtungen. Unter Verwendung dieses Gewichtungsvorgangs konnen Bildkanten nicht-sym- 
metrisch in vier Richtungen verstarkt werden, ohne die vorhergehenden TiefpaBfiltervorgange entlang von 
Kant en sehrim Offset zu storen. Als besonderer Fall, wenn j — ► 4- oo ist, wird nur ein Ausgang von dem Filter mit 
der hochsten Amplitude eff ektiv fur die Kantenverstarkung ausgew&hlt als 

Gleichung (18) 

xh^ni.nz) - xhm(ni^i 2 ), if | Xhro(n t ^i 2 ) | > | Xhn(ni^i 2 ) | 

Es wurde auch festgestellt, dafi die Unterschiedskomponenten di aus der Gleichung (4) anstelle von xw in den 
Gleichungen (16) und (17) verwendet werden konnen, um die Gewichtungsparameter und den endgOltigen 
HochpaBausgang zu bestimmen. Dann ist das hochpaBgefilterte Bild fflr die Kantenverstarkung: 

Gleichung (19) 



Somit ist in diesem Fall viel von der fur das HochpaBfiltern erforderlichen Berechnung eliminiert Es ist jedoch 
anzumerken, da3 nun das HochpaBbild Xbo nicht unabhangig eingestellt werden kann, da di durch die Parameter- 
spezifikationen des TlefpaBfilters bestimmt sind. 

Um die Grenzlinien zwischen unterschiedlichen Bildstrukturen zu verstarken und Oberschwingen von hellen 
Bildkanten in MRI zu verhindern, konnte der Ausgang Xho(ni, 2) von dem obigen richtungsselektiven HochpaBfil- 
ter weiter verarbeitet werden mit einer einfachen SNS (sign non-symmetrical^Funktion, die definiert ist als: 

(Gleichung 20) 

^hotni,n 2 ) = Xho(ni,n2)— w 0 • | xho(tii,n 2 ) | 

wobei w n ein SNS-Gewichtungsparameter ist In einem solchen Fall variieren die SNS-verarbeiteten Hochfre- 
quenzkomponenten xho(n t , n 2 ) mit den Vorzeichen von Xho(ni, 2). Wenn der Parameter w n im Bereich von 0 ^ w n 
< 1 gewahlt ist, dann: 

Gleichung (21) 



Daher kann durch die Wahl von w n eine starkere Kantenverstarkung unter Verwendung eines nicht-symme- 
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trischen HochpaBsignals durchgeftthrt werden, und deutlichere Strukturgrenzlinien kSnnen erhalten werden, 
ohne helle Kanten Uberschwingen einzufQhren. SchlieBlich kann ein Bruchteil dieser SNS-verarbeiteten Fre- 
quenzkomponenten x h0 zu dem tiefpaBverarbeiteten Bild x h o(ni, 2 ) addiert werden, urn das endgflltige rauschge- 
glattete und kantenverstfirkte Bild XhoCm, 2 ) zu erhalten: 6 

5 

Gleichung(22) 

yo(ni,n2) = xio(n lf n 2 )+Wc • Xho(ni,n2) 

to wobei w c der Kantenverstarkungs-Gewichtungsfaktor ist Aus Fig. 4 ist zu sehen, daB dieses Scharfungsschema 
auch eine none RegelmaBigkeit und einen hohen Parallelismus besitzt und daher auch f Qr VLSI- oder Multi-Pro- 
zessorendurchf Qhrung geeignet ist 

Ein Ablauf diagramm fur eine sequentiell durchgefGhrte Computerprogrammdurchf Qhrung des Ausfuhrunes- 
beispiels der vorliegenden Erfindung in einem MRI-System ist in Fig. 5 abgebildet Angenommen, daB die 

is Menuoption "Verfeuierung" gewahlt wurde, dann erfolgt der Einstieg in diese Subroutine bei 500. Auf das 
vorhandene (z. B. herkommliche MR-)Bild wird zugegriffen bei der Subroutine bei 502, und die tiefpaBgefilterten 
UnterbUderm vier Hauptrichtungen werden berechnet bei 504. Danach werden die Unterschiede der Unterbil- 
der im Hmbhck auf das ursprflngliche Filtereingangsbild erzeugt bei 506. Sowohl die Bandbreite des TiefpaBfil- 
terns in Schntt 504 wie auch die bei 508 erzeugten Gewichtungsmuster werden selbstverstandlich durch den 

20 ausgewahlten Satz von Filterparametern (d. h. entweder den fur allgemeine Zwecke ausgewahlten Standardsatz 
Oder einen der neu dehmerten Satze von verfOgbaren Parametern, die vom Bediener als Teil des Auswahimenus 
gewahlt werden konnen) bestimmt Das adaptiv tiefpaBgefilterte Bild wird dann bei 510 synthetisiert Danach 
werden die hochpaBgefilterten Unterbilder in vier Hauptrichtungen bei 512 erzeugt, geeignete Gewichtungsmu- 
ster werden bei 514 erzeugt (wiederum in Abhangigkeit von dem fur die Verwendung ausgewahlten Sate von 

25 Filterparametern). Das adaptiv hochpaBgefilterte Bild wird dann bei 516 erzeugt und mit dem adaptiv tiefpaBge- 
filterten Bild bei 518 kombiniert, urn das endgOItige rauschunterdriickte und strukturverstarkte Filterausganes- 
bild zur Verfugung zu stellen. Danach wird die Subroutine bei 520 verlassen, um zu der MR-Hauptbildverarbei- 
tung zurQckzukehren und die Programm- und MenQwahlmogUchkeiten anzuzeigen. 

Eine gegenwartig bevorzugte AusfOhrung der vorliegenden Erfindung wurde auf einem u-Vax/VMS mit 

30 einem Worner-Matrixprozessor in einem ACCES5™-MRI-System der Fa. Toshiba durchgefflhrt Die Filteroro- 
grammierung ist fur zweidimensionale Bilder geschrieben und kann fur dreidimensionale Bilder erzeugt werden. 
Ein Diagramm dieses Filters ist in Fig. 6 gezeigt 

j w n er L n dem ^nge^nrt 611 F »ter verwendeten Bildmuster haben die gleiche GroBe. Sie haben N Spalten 
und M Reihen. In den meisten FaUen sind N = M0256. Bildmuster sind durch fettgedruckte GroBbuchstaben 
(vne U) bezeichnet, die diesen besonderen Filter im nachfolgenden beschreiben. 

^ i J iS ft M ldes u 311 der SteIle der '" ten SpaIte j " ten Reihe wW be2eichnet durch u <y>' wobei 1 

Bei dem durchgefuhrten Filter werden die Bilder manchmal gefiltert durch eindimensionale Konvolution bzw 
Faltung entlang der vier msbesondere gewahlten Richtungen. Diese vier Richtungen sind die horizontale 
Richtung, die vertikale Richtung, die Hauptdiagonalrichtung und die Nebendiagonalrichtung. Die Filter und die 
geniterten Bildmuster endang dieser Richtungen sind bezeichnet durch die tiefgestellten Buchstaben h. v, d bzw. 
s. Beispielsweise bei einem Kern mit einer Lange 5, etwa 

45 { Q> eo 1^0 a a> \ 

l(1.0+4x©) , (i0+4xa)) , (L0+4x a>f (1.0 +4xfl)} , (L0+4x©)'J 

so definhert durch rCeIIe ^ eindimensionale Faltun e an einem BUd U entlang der horizontalen Richtung 
1.0 



35 



40 



55 



60 



fur jedes ij, und ist bezeichnet durch 
U h « Conv h (U) 



Auch ist zur Vervollstandigung die eindimensionale Faltung an einem BUd U entlang der vertikalen Richtung 
die Hauptdiagonalrichtung und die Nebendiagonalrichtung mit dem vorgegebenen Kern definiert durch 
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u (/./) = 0 ^4 Xfl>) W +-2.y-2)+ ««(/+ w-i)+u(i,j)+ a*i(f-w+i)+ a*J(i-2j +2)) 

fiir jedes i, j und ist jeweils bezeichnet durch 
U v = ConVv(U),Ud = Conva(U), U s - Convs(U) 

Bei der gegenwartigen Durchfuhrung gibt es acht Filterniveaus. Das Filterniveau k kann aus einem Rekon- 
struktionsmenQ gewahlt sein. Refine Jetter: L => P, Q, R, S, T, U f V, W entspricht dem Filterniveau: 1 «= 1, % 3, 4, 5, 
6, 7 bzw. 8. Ein Beispiel fur dieses durchgefQhrte adaptive Filter ist im nachfolgenden fur ein zweidimensionales 
Biid A mit dem Format N mal M und das Filterniveau 1 beschrieben. 

Zuerst werden die vier absoluten lokaJen Unterschiedsmuster durch Bewirken einer eindimensionalen Faltung 
entlang der vier unterschiedlichen Richtungen berechnet, gefolgt vom Erzeugen eines absoluten Wertes fur 
jeden Bildpunkt: 

Dh = Abs/Convh(A)) 
D v ■ Abs/Convh(A)) 
Dd = Abs/Convh(A)) 
D s - Abs/Conv h (A)) 

Der Faltungskern ist der gleiche fiir jede Richtung. Er ist: 



f o ay LO LQ (D o 

l(L0+4xo)) t {L0+4xfl>) t (L0+4xa>)" *(L0+4xa):*(L0+4x4>) 

mit co — 0,5. Der Abs-Operator ist definiert durch: 
V = Abs(U); mit V(U) = |U(iJ)| fur jedes I j. 
AnschlieSend wird ein Schwelienwert r berechnet als definiert durch 



d h. r ist der Durchschnittswert der vier absoluten Unterschiedsmuster. 
Drittens werden die Unterschiedsmuster eingestellt durch: 

D'h mm Max(Dh,a); 
D'v = Max(D v ,a); 
D'd — Max(Dd,a); 
D' s = Max(Ds^T); 

mit 

o = 0,llr, 

wobei 1 das Filterniveau und r die oben berechnete Schweile ist 
Der Max-Operator ist definiert durch: 

V - Max(U,a); mit V(iJ) « max(U(ij)a) fiir jedes i, j. 

Viertens werden die eingestellten Unterschiedsmuster geglattet durch standardmafiige zweidimensionale 
3x3-Faltung: 

D" h - Conv(D' h ); 
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D" v - Conv(D'v); 
D" d - Conv(D'a); 
D" s - Conv(D' s ); 

wobei der Faltungskern in jeder Richtung 



{ co ^LO co 
(h0+2x^)'(L0+2xfl))'(1.0+2XG>) 

undo « 0,25 ist 

(Anmerkung: der vierte Schritt kann weggelassen werden. In diesem Fall ist D"h — D'h, D" v = D' v , D"d «= 
D' dl D" s -D's). 

FQnftens wird ein nicht gerichtetes, gewichtetes hochpaBgefiltertes Bild A" erhalten durch die folgenden 
Schritte: 

a) ein hochpaBgefiltertes Bild A' wird erhalten durch Anwendung einer standardmaBigen zweidimensiona- 
len 5 x 5-Faltung an dem gegebenen Bild A: 

A' - Conv(A), 

wobei der Faltungskern in jeder Richtung 



1 CD CO L0 CO CO 

l(1.0+4x©)'(1.0 + 4x 

undo>= -0,125 x 2/7 = -0,0357 ist. 

b) ein Gewichtungsmuster W wird berechnet durch 

W(U)« 1,0 falls D m (ij)^0^r, 
- Or- D m (ij))x(lr-0,9r) falls 0,9r<D m (U)<lr, 
= 0,0 falls Ir^Dm(ij) 

mit 

Dm(U ) - max(D h (U), Dv(iJ), Dd(ij), Ds(y )) for each i and j. 

c) A" ist der gewichtete Durchschnitt von A und A' unter Verwendung des Gewichtungsmusters W, d. h. 

A"(ij) - A'(ij) + W(ij)[A(iJ).A'(U)]fur jedes i, j. 

(Anmerkung: der f Qnfte Schritt kann weggelassen werden. In diesem Fall ist A"(ij) = A(i j). 
Sechstens werden vier geglattete Bilder aus A" erhalten durch Anwendung einer eindimensionalen Faltung 
entlang der horizontalen, der vertikalen, der hauptdiagonaien bzw. der nebendiagonalen Richtung: 

A"h =< Conv h (A") 
A" v - Convv(A") 
A"d - Conv<j(A") 
A" s - Convs(A") 

Der Faltungskern ist der gleiche for jede Richtung. Er ist 



1 co co jLO c> co 

l(L0 + 4xa>)*(l.O + 4xa)'(L0^ 

mit 

(0 = 0,2+0,1 L 

Siebtens werden vier Gewichtungsmuster 
W" h , W" v , ^" d , W" s berechnet durch 
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W (i i) = W k (i,j) 

41 J wAur+wMr+wAur+wAjy 



*AUJ) +VAi>fr +*AufT +*,iufi* 

Achtens wird ein in vier Richtungen gewichtetes Glattungsbild B als der endgQltige Ausgang berechnet: 
B(iJ) = W h (U)xA''h(ij) + Wv(y)xA'MU)+W d (U)xA-d(U)+Ws(U^ 
furjedesij. 

(Anmerkung: das Ausgangsbild B kann von einem anderen vierfach gerichteten adaptiven HochpaBfilter als 
dem oben beschriebenen verarbeitet werden). 

Einige der in dem gegenwartig durchgef fihrten Filter verwendeten Parameter sind: 

Schritt l:l-d-Faltung mit Kern der Lange 5: 

{1/6,1/6,-4/6,1/6,1/6}; 

Schritt 2: keine 

Schritt3:o = 0,llr;furl - 1.....8; 

Schritt 4: 2-d-Faltung mit Kern der Lange 3 in jeder Richtung: 

{1/6,4/6,1/6} 

Schritt 5: 

a) 2-d-Faltung mit Kern der Lange 5 in jeder Richtung: 
{-0,0416; -0,0416; 1,1666; -0,0416; -0,0416}; 

b) weiche Schwellenfunktion mit unterem cut off Seitenwert: 0,9r und oberem cut off Seitenwert: lr (d. h. ein 
Wert von 1 bis hinauf zu 0,91 und dann ein linearer roll off bis Null bei lr). 

c) keine; 

Schritt 6: 1-d-Faltung mit Kern der Lange 5; 
0,1364; 0,1364; 0,4545; 0,1364; 0,1364} fur! « 1 
0,1538;0,1538;0 r 3846;0,1538;0,1538}fQrl = 2 
0,1667;0,1667;03333;0,1667;0,1667}fQrl - 3 
0,1 765; 0,1 765; 0,2941; 0,1 765; 0,1765} furl « 4 
0,1842;0,1842;0,2632;0,1842;0,1842}furl - 5 
0,1905;0,1905;0,2381;0,1905;0,1905 furl - 6 
0,1957;0,1957;0,2174;0,1957;0,1957 ftirl - 7 
0,2000; 0,2000; 0,2000; 0,2000; 0,2000 furl - 8 
Schritt 7: k - 4 fur Gleichung (8); 
Schritt 8: keine. 

MR-Bilder konnen verwendet werden, urn die Leistung der vorgeschlagenen Rauschglattung, Strukturbewah- 
rung und Filter mit hoher Verarbeitungsgeschwindigkeit zu bewerten. Beispielsweise wurden die Bilder der 
Fig. 7 — 10D erhalten mit einem im Handel erhaltlichen Dauermagnet-lmager (ACCESS™ Toshiba- America 
MRI Inc.) mit geringer Feldstarke (640 G). Alle Bilder in den Experimenten haben die gleiche Bildsignal-Daten- 
anordnungsstruktur 4 {(nl, n2) : 1 nl <> 256, 1 ^ n2 < 256}. 

Aus den Gleichungen (8) und (16) ist zu ersehen, daB groBere Filterparameter k und j die Richtungswahl der 
Filter empfmdlicher fur Bildstrukturen machen, aber groBere Werte fiir k und j erfordern auch mehr Berech- 
nung. Auch gibt es einen Abgleich zwischen Kantenbewahrung und Rauschglattung, wenn der TiefpaBparame- 
ter k gewahlt wird. Wenn k grSBer wird, wird die Kantenbewahrungsleistung des Filters besser, wahrend die 
RauschglSLttungsleistung schwacher wird. In der Praxis stellt ein Filter mit einem Wert k « 4 gute Resultate fQr 
Rauschglattung und Kantenbewahrung zur Verftigung. 

Bei dem ersten Experiment wurde eine sagittale Kopfabbildung aus einer dreidimensionalen (3D)-Erfassiing 
verwendet. Das ursprungliche Bild ist in Fig. 7 gezeigt Vier ausgerichtete TiefpaBfilter, wie in Gleichung (2) 
definiert mit ai = 1/5 wurden verwendet, urn die vier in Fig. 8A, 8B, 8C bzw. 8d gezeigten Unterbilder zu 
erhalten. Der exponentietle Parameter k — 4 wurde fiir die Bestimmung der Gewichtungsparameter in Glei- 
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chung (8) verwendet. Diese Unterbilder haben starke ausgerichtete Strukturen als Ergebnis der ausgerichteten " 
TiefpaBfilterung. Kanten- und Linienstrukturen entlang einer Richtung sind geglattet und bewahrt durch das 
Filter in der gleichen Richtung, aber verschwommen durch die Filter in anderen Richtungen. Die Komponenten 
jedes Teilbildes sind in Abhangigkeit von den Bildunterschieden gewahlt, um zu dem endgiiltigen tiefpaBgefilter- 

5 ten und kantenverstarkten Bild, das in Fig. 9 gezeigt ist, beizutragen. Es ist klar ersichtlich, daB Rauschen im BHd 
der Fig. 7 effektiv reduziert worden ist und daB die Bildstrukturen in Fig. 9 gut bewahrt worden sind. Rauschre- 
duzierung kann auch verstanden werden durch den Vergleich von Rausch-SD-Werten vor und nach dem Filtern. 
Der SD-Wert wurde geschatzt an einem nicht-texturierten Biidbereich mit 200 Bildpunkten (z. B. siehe L. 
Kaufmann, D.M. Kramer, L.E. Crooks und DA. Ortendahl, "Measruring Signal-to-Noise-Rations in MR Ima- 

10 ging", Radiology 1989; 173: 265—67). Der SD-Wert in Fig. 7 betragt 83 und ist in Fig. 9 durch Filtern auf 42 
reduziert Zwei transaxial Kopfabbildungen (eine basierend auf den ersten auftretenden NMR-Spinechoreak- 
tionen und eine basierend auf den zweiten auftretenden NMR-Spinechoreaktionen) aus einer zweidimensiona- 
len (2D) Erfassung wurden auch verwendet, um die Leistungsfahigkeit des Filters besser verstehen zu konnen. 
Diese zwei Bilder haben geringere Rauschabstande als das vorherige 3D-Bild, insbesondere das Zweitechobild. 

15 Das ursprttngliche Bild und das gefilterte Bild fur das erste Echo sind in Fig. 10A und Fig. 10B gezeigt. Das 
ursprungliche Bild und das gefilterte Bild fur das zweite Echo sind in Fig. IOC bzw. Fig. 10D gezeigt Wieder 
wurden 200 Bildpunkte fur SD-Schatzungen verwendet Es ist zu sehen, daB die gefilterten Bilder eine viel 
bessere visuelle Qualitat als die ursprunglichen Bilder haben. 

Basierend auf der beschriebenen Filterstruktur wurde ein Matrbcprozessor (Warrior, SKY Computers, Inc.) 

20 verwendet, um diesen Algorithmus in Verbindung mit einem Micro VAX-I-System (Digital Equipment Corpora- 
tion) durch geeignete Programmierung eines ACCESS™-MRI-Systems von Toshiba (das schon solche Hardwa- 
rekomponenten beinhaltet) durchzufuhren. Eine Synchronanordnung von parallelen Verarbeitungselementen 
(PE, processing elements) im Matrizenrechner bewaltigt einen Befehl und mehrfache Daten, um eine hone 
Verarbeitungsgeschwindigkeit zu erreichen. Als Ergebnis wird ein hoher Durchsatz bei der Bildverarbeitung 

25 erreicht Die Zeit, die zur Verarbeitung eines Bildes mit 256 mal 256 Bildpunkten benotigt wird, kann weniger als 
5 Sekunden betragen. Da vier Unterbilder in diesem Filter gespeichert werden mussen, ist ein Nachteil dieses 
Filters im Vergleich mit koriventionellen adaptiven Filtern, die lokale Operationen oder rekursive Ansatze 
verwenden. die Notwendigkeit, mehr Bildspeicherraum zur Verfugung zu stellen. Dieses Filter ist auch geeignet 
fur die Durchfahrung mit Mehrrechner-Systemen. Anders als Matrizenrechner beinhaltet ein Mehrrechner-Sy- 

30 stem mehrere Prozessoren, und diese Prozessoren kommunizieren und kooperieren auf unterschiedlichen 
Niveaus zur Losung eines gestellten Problems (z. B. siehe K. Hwang, Computer Architecture and Parallel 
Processing, McGraw-Hill, Ino, 1984, auf das hiermit vollinhaltlich Bezug genommen wird). Aus Fig. 3 und 4 ist zu 
sehen, daB sowohl der TiefpaBfilterteil und der Kantenscharf ungsteil dieses Filters hochparallel sind. Daher kann 
die Berechnungsaufgabe, die bei vier unabhangigen Hoch-/TiefpaBfilterkanalen anfallt, aufgeteilt und verschie- 

35 denen Prozessoren zugeteilt werden. Wenn beispielsweise das Mehrrechner-Computersystem von Mercury 
Computer System (z. B. siehe J. Goddard, J. Greene, B.S. Isenstein, F.P. Lauginiger und L.C. Walsh, "Multicompu- 
ter Implementations of Reconstruction and Postprocessing Algorithms , ^ Proceedings of SMRM 11th Annual 
Meeting, Berlin, August 1992) verwendet wird, kann es jedem Prozessor in diesem Vierprozessor-System 
zugeteilt werden, xu, dj, wu und xu auf der TiefpaBstuf e und whu whj und hi auf der HochpaBstuf e zu berechnen, um 

40 eine parallele Verarbeitung fQr schnelle Bildverstarkung zu erreichen. Die Datenkommunikation zwischen 
unterschiedlichen Kanaien ist nur erforderlich beim Berechnen der Gewichtungsparameter w'us und w'ws, 
weshalb eine Hochgeschwindigkeitsverarbeitungzu erwarten ist 

Ein einfaches, schnelles Bildfilter wurde vorgestellt Das Filter reduziert effektiv Bildrauschen, wahrend es 
Bildstrukturen bewahrt Die Filterleistung und die Berechnungseffizienz machen es nutzlich fur schnelle Rau- 

45 schreduzierung von medizinischen Bildern und damit auch hilfreich fur viele weitere Korrekturen von Bildver- 
zerrungen. Die parallele Struktur des Filters kann verwendet werden, um eine weiter erhohte Verarbeitungsge- 
schwindigkeit durch die Verwendung von Multiprozessor-Systemkonfigurationen zu erreichen. 

Patentanspriiche 

50 

1. Verfahren zum Filtern von eingegebenen digitalen Bildsignaien, welche die angeordneten Pbcel eines 
digitalisierten Bildes darstellen, wobei das Verfahren die Schritte aufweist: 

Erzeugen einer Mehrzahl von digitalisierten Satzen von Unterbildsignalen durch TiefpaBfiltern der einge- 
gebenen digitalen Biidsignale entlang einer jeweiligen entsprechenden Mehrzahl von Richtungen innerhalb 
55 des Bildes; und 

Erzeugen eines Satzes von rauschgeglatteten digitalen Bildsignalen als eine gewichtete Kombination der 
Mehrzahl von Satzen von Unterbildsignalen. 

2. Verfahren nach Anspruch 1, wobei die gewichtete Kombination jedes Satzes von Unterbildsignalen 
durchgef uhrt wird in Abhangigkeit des Unterschiedes zwischen dem Unterbild und dem eingegebenen Satz 

60 von digitalen Bildsignalen. 

3. Verfahren nach Anspruch 2, wobei der Beitrag eines Unterbiides zu dem rauschgeglatteten Bild mit 
zunehmendem Unterschied abnimmt 

4. Verfahren nach Anspruch 1 , welches des weiteren aufweist: 

Erzeugen einer Mehrzahl von weiteren Satzen von digitalisierten Unterbildsignalen durch HochpaBfiltern 
65 des rauschgeglatteten digitalen Bildsignals entlang einer jeweils entsprechenden Mehrzahl von Richtungen 

innerhalb des Bildes; und 

Erzeugen eines Satzes von gescharften digitalen Bildsignalen als eine gewichtete Kombination der Mehr- 
zahl von weiteren Satzen von Unterbildsignalen. 
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5. Verfahren nach Anspruch 2, wobei eine gewichtete Kotnbination der Unterschiede zwischen den Unter- 
bildern und dem eingegebenen digitalen Bildsatz verwendet wird, urn einen Satz von gescharften digitalen 
Bildsignalen herzustellen. 

6. Verfahren nach Anspmch 4 oder 5, wobei die rauschgegeglatteten digitalen Bildsignale und ein Abschnitt 

der gescharften digitalen Bildsignale kombiniert werden, urn einen zusammengesetzten gefilterten Satz von 5 
Ausgangsbildsignalen herzustellen. 

7. Verfahren nach Anspruch 1 oder 4, wobei die TiefpaBfilterung durch LSi-(linear space invariant)-TiefpaB- 
filtem durchgefuhrt wird 

8. Verfahren nach Anspruch 7, bei dem das LSi-TiefpaBfiltern in jeder der Richtungen durchgefuhrt wird 
durch Erzeugen fur jedes Pixel einer gewichteten Kombination einer vorbestimmten Anzahl seiner nach- 10 
sten Nachbarn in einer Linie parallel zu dieser Richtung. 

9. Verfahren nach Anspruch 1 oder 4, bei dem die Mehrzahl von Richtungen des TiefpaBfilterns vier 
Richtungen umfaBt, wetche im wesentlichen in 45° - Winkelintervallen innerhalb des Bildes angeordnet sind 

10. Verfahren nach Anspruch 8, bei dem die Mehrzahl von Richtungen des TiefpaBfilterns vier Richtungen 
umfaBt, welche im wesentlichen in 45° - Winkelintervallen innerhalb des Bildes angeordnet sind 15 

11. Verfahren nach Anspruch 4, bei dem das TiefpaBfiltern durch nicht-syrnmetrisches HochpaBfiltern 
durchgefuhrt wird 

12. Verfahren nach Anspruch 1 1, bei dem nicht-symmetrisches HochpaBfiltern in jeder der Richtungen fur 
HochpaBfiltern durchgefQhrt wird durch Erzeugen fur jedes Pixel einer gewichteten Kombination einer 
vorherbestimmten Anzahl seiner nachsten Nachbarn in einer Linie parallel zu dieser Richtung. 20 

13. Verfahren nach Anspruch 12, bei dem die gewichtete Kombination nicht-symmetrisch im Hinblick 
sowohl auf die raumliche Richtung als auch auf die Polaritat erzeugt wird 

14. Verfahren nach Anspruch 12 oder 13, bei dem die Mehrzahl von Richtungen des HochpaBfil terns vier 
Richtungen umfaBt, welche im wesentlichen in 45° -Winkelintervallen innerhalb des Bildes angeordnet sind 

15. Vorrichtung zum Filtern von. eingegebenen digitalen Bildsignalen, welche die angeordneten Pixel einer 25 
digitalen Abbiidung darstellen, wobei die Vorrichtung aufweist: 

eine Einrichtung zum Erzeugen einer Mehrzahl von Satzen von digitalisierten Unterbildsignalen durch 
TiefpaBfiltern des eingegebenen digitalen Bildsignais entlang einer jeweils entsprechenden Mehrzahl von 
Richtungen innerhalb des Bildes; und 

eine Einrichtung zum Erzeugen eines Satzes von rauschgeglatteten digitalen Bildsignalen als gewichtete 30 
Kombinadon der Mehrzahl von Unterbildsignalsatzen. 

16. Vorrichtung nach Anspruch 15, welche eine Einrichtung aufweist, die verursacht daB die gewichtete 
Kombination jedes Unterbildsignalsatzes in Abhangigkeit des Unterschiedes zwischen dem Unterbild und 
dem eingegebenen digitalen Bildsignalsatz aixsgeffihrt wird 

17. Vorrichtung nach Anspruch 16, welche eine Einrichtung aufweist, die verursacht, daB der Beitrag eines 35 
Unterbildes zu der rauschgeglatteten digitalen Abbiidung abnimmt, wenn der Unterschied zunimmt 

1 8. Vorrichtung nach Anspruch 1 5, welche des weiteren aufweist: 

eine Einrichtung zum Erzeugen einer Mehrzahl von weiteren Satzen digitalisierter Unterbildsignale durch 
HochpaBfiltern des rauschgeglatteten digitalen Bildsignais entlang einer jeweils entsprechenden Mehrzahl 
von Richtungen innerhalb der Abbiidung; und 40 
eine Einrichtung zum Erzeugen eines Satzes von gescharften digitalen Bildsignalen als eine gewichtete 
Kombination der Mehrzahl weiterer Satze von Unterbildsignalen* 

19. Vorrichtung nach Anspruch 16, welche, eine Einrichtung aufweist, die bewirkt daB eine gewichtete 
Kombination der Unterschiede zwischen den Unterbildern und dem eingegebenen Satz von digitalen 
Bildern zum Herstellen eines Satzes von gescharften digitalen Bildsignalen verwendet wird 45 

20. Vorrichtung nach Anspruch 18 oder 19, welche eine Einrichtung aufweist, die bewirkt, daB die rauschge- 
glatteten digitalen Bildsignale und ein Abschnitt der gescharften digitalen Bildsignalen zur Herstellung 
eines zusammengesetzten, gefilterten Satzes von ausgegebenen Bildsignalen kombiniert werden. 

21. Vorrichtung nach Anspruch 15 oder 18, welche eine Einrichtung zuin Bewirken des TiefpaBfilterns durch 
LSI-HefpaBfiltern aufweist 50 

22. Vorrichtung nach Anspruch 21, welche eine Einrichtung zum DurchfUhren des LSI-HefpaBfilterns in 
jeder der Richtungen durch Erzeugen fur jedes Pixel einer gewichteten Kombination einer bestimmten 
Anzahl ihrer nachsten Nachbarn auf einer Linie parallel zu dieser Richtung aufweist 

23. Vorrichtung nach Anspruch 15 oder 18, welche eine Einrichtung aufweist, die bewirkt, daB die Mehrzahl 
von Richtungen des TiefpaBfilterns vier Richtungen umfaBt, die in Winkelabstanden von im wesentlichen 55 
45° innerhalb des Bildes ausgerichtet sind 

24. Vorrichtung nach Anspruch 22, welche eine Einrichtung aufweist die bewirkt daB die Mehrzahl von 
Richtungen des TiefpaBfilterns vier Richtungen umfaBt, die in Winkelabstanden von im wesentlichen 45° 
innerhalb des Bildes ausgerichtet sind 

25. Vorrichtung nach Anspruch 18, welche eine Einrichtung zum Durchftihren des HochpaBfil terns durch 60 
nicht-symmetrisches HochpaBfiltern aufweist 

26. Vorrichtung nach Anspruch 25, welche eine Einrichtung zum Durchf Qhren des mcht-symmetrischen 
HochpaBfUtems in jeder der Richtungen aufweist durch Erzeugen fur jedes Pixel einer gewichteten Kombi- 
nation eirier bestimmten Anzahl seiner nachsten Nachbarn auf einer Linie parallel zu dieser Richtung. 

27. Vorrichtung nach Anspruch 26, welche eine Einrichtung zum nicht-symmetrischen Erzeugen der gewich- 55 
teten Kombination im Hinblick sowohl auf die raumliche Richtung als auch auf die Polaritat aufweist. 

28. Vorrichtung nach Anspruch 26 oder 27, welche eine Einrichtung aufweist, um zu bewirken, daB die 
Mehrzahl von HochpaBfilterrichtungen vier Richtungen umfaBt, die in Winkelabstanden von im wesentii- 
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chen 45° innerhalb des Bildes ausgerichitet sind. 
29 * MRI-System mit einem MR-Abbilduhgsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet, 
wobei der MR-Abbildungsdaten-Prozessor so programmiert ist, daB er das FUterverfahren von Anspruch 1 
bewirkt r 



30. MRI-System mit einem MR-AbbUdungsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet, 
wobei der MR-Abbildungsdaten-Prozessor so programmiert ist, daB er das FUterverfahren von Anspruch 2 
bewirkt 

31 • ^M-System mit einem MR-Abbildungsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet, 
bewirkt MR " Abbildungsdaten " Prozessor so Programmiert ist, daB er das Filterverfahren von Anspruch 3 

32. MRI-System mit einem MR-Abbildungsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet, 
wobei der MR-Abbildungsdaten-Prozessor so programmiert ist, daB er das Filterverfahren von Anspruch 4 
bewirkt 

33. MRI-System mit einem MR-Abbildungsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet, 
wobei der MR-Abbildungsdaten-Prozessor so programmiert ist, daB er das Filterverfahren von Anspruch 5 
bewirkt K 

34. MRI-System mit einem MR-Abbildungsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet 
wobei der MR-Abbildungsdaten-Prozessor so programmiert ist, daB er das Filterverfahren von Anspruch 6 
bewirkt 

35. MRI-System mit einem MR-Abbildungsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet, 
wobei der MR-Abbildungsdaten-Prozessor so programmiert ist, daB er das FUterverfahren von Anspruch 7 
bewirkt 

36. MRI-System mit einem MR-Abbildungsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet, 
bewSct R " Abbl,dungsdaten " Prozessor so Programmiert ist, daB er das FUterverfahren von Anspruch 8 

37. MRI-System mit einem MR-Abbildungsdaten-Prozessor, welcher einen Matrbcprozessor beinhaltet, 
wobei der MR-Abbildungsdaten-Prozessor so programmiert ist, daB er das FUterverfahren von Anspruch 9 



HierzulO Seite(n) Zeichnungen 



16 



Nummer: DE 44 31 343 A1 

Int. CI.*: H04N 5/213 

Offeniegungstag: 9. Marz 1995 




408070/644 



ZE1CVWUNGEN SEVTE 2 



Nummer: 
Int. CI.*: 

Offentegungstag: 



OB 44 31 349 At 
H04N 5/213 

9.Marz1995 * 



8 




CM 
CO 



to 



u 




o 




w 




03 




O 


U 


S 


o 




1 0) 




to 


a 


•H (U 


1 o 


M N 


H rH 


•p o 


S CQ 


X 04 




O 
CO 



00 
-CM 



N 

I 

O 

04 

p 

<H 
*H 



C 

O 
P 

u 

O 
-P 



O 
CP 
! 

>H C 
«-J O 



C 

o x: 

> O 
03 -Q 

O M <p 

cr»^: p 
S3 o c 
o 6 :=> 

N 



Q 

s 
o 

v; 

c c 

P CO 

-c 3 
o 03 
*H M 
£ 03 

<D *CJ 

C 
^ O 
O *H 
C P 
•H 03 
<U C 



-C c 

O 0) 

TO -p 
•H 

V* <D 

-P 4J 

0 rH 

1 *H 
>i <D 

to 

■P C r-4 

-C O H 

o aj a 

•H «W TO 
C M 

:o3 _ 

6 .C 03 

GJ U ^ 

C CO TO 

'HOD 



W 03 



408 070/644 



ZEICHNUNGEN SEITE 3 



Nummer: 
Int. CI. 6 : 

Offenlegungstag: 



DE 44 31 349 A1 
H04N 5/213 

9. Ma rz 1995 



90° 




FIG.2 



408 070/644 



ZBCHNUNGEN SEITE 4 



Nummer; 
Int. CI. 6 : 

Offenlegungstag: 



DE 4431 349 AJ 
H04N 5/213 

9.Marz1995 * 



FIG. 4 



*io ^ 



I 


i 


I 


1 


fhi(xto) 




Aa(*io) 












{ 



408 070/644 



ZCICHNUNGEN SEITE 5 M 

Nummer: DE44 31349 A1 

lnt. CI.6: H 04N 5/213 

Offenlegungstag: 9. Marz1995 



c 



I 



I 



I 



1 



I 



I 



I 



I 



I 



3 

3 



500 



502 

Beginnen mit dem vorgegebenen 
x 0 
504 



] 



FIG. 5 

Ber. der Unterbilder x 
durch TiefpaBf iltern vo£ 
Xq in vier Hauptrichtungen 



J 



506 

Berechnen der Uhterschiede d 
zwischen den Unterbildern x 1 
xma dem urspruglichen Bila x A 



r 508 



Berechnen der Gesichtungsmuster 



I 



510 



Synthetisieren des adaptiv 
tiefpaBgefilterten Bildes 



10 
512 



I Ber. der Unterbilder x 
BachpaSf iltern von x. J"- 
in vier Hauptrichtun 15 



durch 
in 



Ber. der Gewichtungsmuster w h± 



r516 
Syn 
hoc 



Synthetisieren des adaptiv 
hochpaGgef i 1 terten Bildes 5^ 



] 



3 



r sis 

Ersrallen des enagiilrig 
rauschunrerdrucKten und 
strukturverstarkten 
Bildes y durch Kbmbiniere: 



520 



von x 1Q und ^ 



Ausgang zu Haupt-ivjRi-Bildverarbeit- 
ungs una Anzeigeprograimen 



408 070/644 



ZE1CHNUNGEN SEITE 6 



Nurnmer: 
Int. CI. 6 : 

Offenlegungstag: 



DE 4431 349 A1 
H04N 5/213 
9. Marz1995 * 




FIG. 6 



i, D v °d 



D 'h D 'v D 'd P's | 



D"h D" v D"d D" 



s 



8 



a 



Wh W Wd w s 



I 



408 070/644 



2EICHNUNGEN SEtTE 7 



Nummer: 
Int. CI. 6 : 

Offenlegungstag: 



DE 44 31 349 A1 
H04N 5/213 

9. Marz1995 




408070/644 



ZttCHNUNGEN SETTE a Nummer: D£ 4431 349 A1 

Int. CI.*: H 04 N 5/213 

Offenlegungstag: 9. Merz1995 * 




FIG. 8c 



FIG. 8d 




408070/644 



ZEICHNUNGEN SEITE 9 



Nummer: 
Int. CI. 6 : 

Offenlegungstag: 



DE 4431 349 A1 
H04N 5/213 
9. Marz1995 




FIG. 9 



408070/644 



ZEICHNUNGEN SEITE 10 



Niimroer: 
Int. CI. 6 : 

Off enleg ungstag : 



OE 44 31 349 AT 
H04N 5/213 

9.Marz1995 # 



FIG. 10a 




FIG. 10c 





-.7? 




._. .. 1 . * 


• * 

■ \ 





FIG. 10b 




408070/644 



This Page is Inserted by IFW Indexing and Scanning 
Operations and is not part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of the original 
documents submitted by the applicant. 

Defects in the images include but are not limited to the items checked: 

□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SIDES 
i^FADED TEXT OR DRAWING 
^BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 
tTREFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 



□ OTHER: 



IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 



